Methods and compositions for tagging and analyzing samples

ABSTRACT

The invention relates to methods of tagging analytes in a sample.

CROSS-REFERENCE

This application claims the benefit of U.S. Provisional Application No. 61/801,785, filed Mar. 15, 2013, and U.S. Provisional Application No. 61/806,143, filed Mar. 28, 2013, both of which applications are incorporated herein by reference in their entireties.

BACKGROUND OF THE INVENTION

The analysis of nucleic acids and proteins in biological samples is an essential element of molecular biology. The ability to detect, discriminate, and utilize genetic and proteomic information allows sensitive and specific diagnostics as well as treatment.

The present invention provides for rapid tagging and analysis of nucleic acids and protein at the single cell level.

SUMMARY OF THE INVENTION

The present invention provides methods of tagging target oligonucleotides that include steps of (a) partitioning DNA into a plurality of compartments; (b) performing an in vitro transcription reaction on the DNA within the compartments, thereby obtaining compartments comprising RNA; (c) merging the interior of the compartments comprising RNA with the interior of a set of compartments comprising target oligonucleotides; (d) hybridizing the RNA to the target oligonucleotides; and (e) performing a reaction to attach a sequence corresponding to the RNA to the target oligonucleotides.

In some embodiments, the DNA is double stranded.

In some embodiments, the compartments are droplets within an oil-and-water emulsion.

In some embodiments, the target oligonucleotides include at least one target oligonucleotide comprising a cell tag and a molecule tag.

In some embodiments, the target oligonucleotides are DNA.

In some embodiments, the method may further include a step of partitioning the target oligonucleotides into the set of compartments before the merging of step (c).

In some embodiments, the method may further include a step of partitioning a set of cells into the set of compartments and lysing the cells in order to liberate cellular oligonucleotides before the merging of step (c). For example, the cellular oligonucleotides are the target oligonucleotides. Alternatively, the cellular oligonucleotides are cellular mRNA.

In some embodiments, the method may further include a step of conducting a reverse transcription on the cellular mRNA to generate cellular cDNA.

For example, the reverse transcription reaction is performed with a primer specific for a region of the genome. The region of the genome can be an immunoglobulin gene or a T-cell receptor gene.

For example, the reverse transcription reaction can be conducted within the set of compartments prior to the merging step.

Alternatively, the reverse transcription reaction can be conducted within the merged compartments.

In some embodiments, the target oligonucleotides are cellular cDNA.

In some embodiments, the reaction is a Rapid Amplification of cDNA Ends (RACE) reaction.

In some embodiments, the DNA is conjugated to a solid support. For example, the solid support is a bead.

The present invention also provides methods that include steps of (a) providing a plurality of beads comprising a plurality DNA oligonucleotides; (b) providing a plurality DNA oligonucleotides comprising a primer sequence, a universal sequence an adapter sequence and a cellular tag; (c) merging the beads of step (a) and the oligonucleotides of step (b) into a plurality of compartments such that each compartment comprises a single bead and a single oligonucleotide; (d) performing an amplification reaction on the oligonucleotides within the compartments, thereby obtaining a plurality of DNA oligonucleotides comprising the primer sequence, the universal sequence, the adapter sequence and the cellular tag; (e) performing an in vitro transcription reaction on the DNA within the compartments, thereby obtaining compartments comprising RNA the primer sequence, the universal sequence and cellular tag; (f) merging the interior of the compartments comprising RNA with the interior of a set of compartments comprising target oligonucleotides; (g) hybridizing the RNA to the target oligonucleotides; and (h) performing a reaction to attach a sequence corresponding to the RNA to the target oligonucleotides.

In some embodiments, the plurality of oligonucleotides on the bead includes a molecule tag.

The present invention further provides methods of tagging target oligonucleotides that include steps of (a) isolating a plurality of mRNA from a biological sample comprising a plurality of cell types; and (b) performing reverse transcription of the mRNA using a primer specific for the target oligonucleotide and a template switching oligonucleotide comprising a molecule tag, a universal sequence, and an adapter sequence to produce tagged target cDNA.

In some embodiments, the target cDNA is tagged at the 3′end.

In some embodiments, the target oligonucleotide is an immunoglobulin or T-cell receptor.

In some embodiments, the adapter sequence is specific to a sequencing platform.

In some embodiments, the molecule tag is an oligomer. For example, the oligomer is a randomer.

In some embodiments, the randomer is at least a 9 mer.

In some embodiments, the method further includes a step of amplifying the target cDNA using the universal sequence and a primer specific for the target oligonucleotide.

In some embodiments, the method may further include a step of sequencing the amplified cDNA.

Also provided are methods of determining the immune repertoire in a subject by (a) isolating a plurality of mRNA from a biological sample comprising a plurality of cell types; (b) performing reverse transcription of the mRNA using a immunoglobulin or T-cell receptor specific primers and a template switching oligonucleotide comprising a molecule tag, a universal sequence and an adapter sequence to produce molecule tagged immunoglobulin or T-cell receptor cDNA; (c) amplifying the cDNA using the universal sequence and a primer specific for the target oligonucleotide; (d) sequencing the cDNA to produce a plurality of sequencing reads; (e) grouping the sequence reads with the same molecule tag and clustering the sequences within the same group; and (f) building a consensus sequence for each cluster to produce a collection of consensus sequences wherein the consensus sequence is used to determine the diversity of the immune repertoire.

In some embodiments, the target cDNA is tagged at the 3′end.

INCORPORATION BY REFERENCE

All publications, patents, and patent applications mentioned in this specification are herein incorporated by reference to the same extent as if each individual publication, patent, or patent application was specifically and individually indicated to be incorporated by reference.

BRIEF DESCRIPTION OF THE DRAWINGS

The novel features of the invention are set forth with particularity in the appended claims. A better understanding of the features and advantages of the present invention will be obtained by reference to the following detailed description that sets forth illustrative embodiments, in which the principles of the invention are utilized, and the accompanying drawings of which:

FIG. 1 is a schematic overview of the performance of an assay on a patient sample.

FIG. 2 is a schematic overview of an assay for tagging cells within a sample.

FIG. 3 depicts steps that may be involved in a Rapid Amplification of cDNA Ends (RACE) assay.

FIG. 4 depicts additional steps that may be involved in a RACE assay.

FIG. 5 depicts additional steps that may be involved in a RACE assay.

FIG. 6 depicts performance of an immune repertoire assay on a population of cells.

FIG. 7 depicts flow cytometry and gene expression assays for the immune repertoire.

FIG. 8 depicts tagging and analysis of analytes within partitions that may be merged.

FIG. 9 depicts multiplexed tagging and analysis.

FIG. 10 depicts preparation of an oligonucleotide-conjugated bead and generation of a tag-bearing oligonucleotide from the bead.

FIG. 11 depicts preparation of an oligonucleotide-conjugated bead and generation of a tag-bearing oligonucleotide from the bead.

FIG. 12 depicts preparation of an oligonucleotide-conjugated bead and generation of a tag-bearing oligonucleotide from the bead.

FIG. 13 depicts preparation of an oligonucleotide-conjugated bead and generation of a tag-bearing oligonucleotide from the bead.

FIG. 14 shows steps for tagging and analyzing mRNA.

FIG. 15 shows additional steps for tagging and analyzing mRNA.

FIG. 16 depicts examples of methods of partitioning and merging.

FIG. 17 depicts combinatorial properties for certain randomer tag lengths.

FIG. 18 depicts the use of information systems in analysis.

DETAILED DESCRIPTION OF THE INVENTION I. Overview

This disclosure provides methods and compositions for tagging molecules and subsequent analysis of the tagged molecules. In some cases, the disclosure provides methods and compositions for partitioning analytes (e.g., cells, polynucleotides) into individual partitions (e.g., droplets, wells, spots on arrays, etc.) and also provides methods and compositions for labeling the analytes within the partition.

In some cases, this disclosure generally relates to obtaining a sample (e.g. blood, saliva, tissue, cells) from a subject [110] (e.g. human, animal, plant, fungus, bacteria, population of cells, biofilm), partitioning and labeling the sample components with tagging molecules [120], processing and analyzing the labeled sample components [130], and/or reporting results from the analysis [140]. Cellular components such as DNA and RNA, individual cells or a population of cells may be partitioned. Cell population can include cells of similar function, such as for example, immune cells (e.g. B-cells or T-cells), cancer cells, or nerve cells. Moreover, a population of cells may be partitioned into cell partitions [210] and, separately, a population of beads with surface-bound oligonucleotides (e.g., DNA, dsDNA) may be partitioned within bead partitions [230]. The individual cells may be lysed within the cell partitions [220]. In some cases, the oligonucleotides are transcribed from the beads [240] such that RNA transcripts are produced in the bead partitions. In some cases, the RNA transcripts are reverse-transcribed into cDNA within the bead partitions; in some cases, the RNA transcripts are reverse-transcribed into cDNA within the bead partitions at a later step, such as after the contents of the bead partitions (e.g., droplets) are merged with the contents of another partition (e.g., a different droplet). The cell partitions may be combined with the bead partitions on a partition-by-partition basis [250]. Tagging reactions may then be conducted within the combined partitions in order to tag the cellular mRNA (or cDNA derived therefrom) with the oligonucleotides derived from the beads [250]. The tagged products may then be pooled, amplified, and sequenced [260].

The methods and compositions provided in this disclosure allow for tracking information sources and preserving the heterogeneity of information in a sample as it is analyzed. By labeling a sample at the individual component level, the resolution of the information can be maintained at the individual component level throughout the stages of analysis, regardless of subsequent merging or combining of the components. For example, a sample containing multiple cell types can be partitioned to a single cell type. Cell-types that can be partitioned into a single cell type include for example immune cells such as B-cells or T-cells. Alternatively, a sample containing multiple cell types can be partitioned into individual partitions containing single cells. By labeling the information-bearing molecules (e.g. DNA, RNA, protein) from a single cell type or a single cell, the individual partitions can then be merged for further analysis without loss of the single cell-level information resolution. The methods and compositions provided herein may also include additional labels, such as labels that enable quantification of the individual molecules within a partition. For example, a single partition may contain numerous unique labels, each with a different sequence that can be used to tag and quantify individual molecules within a partition.

The methods, compositions and kits provided herein are broadly applicable to a variety of life science-related fields, including biomedical research, drug discovery and development, and clinical diagnostics. Potential applications include gene expression profiling at the single cell type or single cell level for the detection and/or monitoring of cancer, autoimmune disease, viral infection, organ transplant rejection, and other diseases or disorders. The present disclosure may also be used to (a) analyze the immune repertoire of a subject, such as a subject with a particular disease or disorder; (b) elucidate intracellular signaling pathways; (c) validate therapeutic targets for drug discovery and development; and (d) identify or detect biomarkers, particularly biomarkers related to normal or diseased biological states. The present disclosure may also be used to analyze circulating cell-free DNA or RNA in order to predict, monitor, detect and/or diagnose conditions or diseases, including organ rejection.

The methods and compositions disclosed herein offer several important advantages over existing techniques for monitoring gene expression in cells or tissue. Importantly, the methods and compositions provided herein enable the monitoring and detection of gene expression in single cells, thereby eliminating the systematic errors and noise that may arise due to sampling of heterogeneous cell populations when collecting data using conventional techniques. The existence of heterogeneous cell populations in test samples may arise, for example, from asynchronous cell division in populations of cultured cells; or, in some cases, heterogeneity may be due to mixtures of different cell types present in tissue samples, biofilms, bioreactor samples, blood samples, biopsy samples or other complex samples. Another important advantage of the methods and compositions disclosed herein is the potential for eliminating or reducing errors caused by PCR amplification bias, for example, through the use of molecular tags that label different molecules within a sample. For example, if, when analyzed, the same molecular sequence is found to have two different tags, this may indicate that there were two copies of the molecule within the partition. This information may also be useful to discount results due to amplification errors. A third important advantage of the methods and compositions disclosed herein is the potential for expanding the range of biomarkers used to sort and classify cells. In addition to targeting gene sequences that code for the extracellular protein markers, the approach described herein enables the use of intracellular markers, for example gene sequences coding for transcription factors or cytokines, for cell sorting and classification in order to facilitate correlations between gene expression and cell function. The present disclosure thus offers an approach for obtaining higher quality genomic data from biological samples, and thus the potential for developing better therapeutics and improved detection of disease.

II. Assays

A. Labeling within Partitions

This disclosure provides methods and compositions for tagging analytes at a single component level, such as at the level of a single cell type or single cell. In some cases, analytes are partitioned into a set of partitions; labels are partitioned into a separate set of partitions, and the contents of individual partitions within each set are combined to enable labeling of the analyte.

Partitioning Analytes

A sample comprising analytes (e.g., cells) can be partitioned into a set of individual partitions (e.g., droplets or wells). In some cases, a partition within the set of individual partitions contains at most one analyte. In some cases, a partition within the set of individual partitions contains at most 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, 30, 40, 50, 60, 70, 80, 90, or 100 analytes. In some cases, a partition within the set of individual partitions contains, on average, one analyte. In some cases, a partition within the set of individual partitions contains, on average, 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, 30, 40, 50, 60, 70, 80, 90, or 100 analytes, or more. In some cases, the set of individual partitions comprises empty partitions. Often, the set of individual partitions comprises some empty partitions and some partitions comprising analytes (e.g., at most one analyte, at most two analytes, etc.). In some cases, an analyte comprises a plurality of components (e.g., a plurality of molecules). In some cases, a technique is applied to ensure that all of the partitions comprise at most one analyte; for example, the empty partitions (e.g., droplets) may be sorted out by a flow sorter.

Partitioning Tags

Tags or solid supports (e.g., beads) conjugated to tags can be partitioned into a set of individual partitions. In some cases, a partition within the set of individual partitions contains at most one tag or solid support (e.g., a bead) conjugated to a tag. In some cases, a partition within the set of individual partitions contains at most 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, 30, 40, 50, 60, 70, 80, 90, or 100 tags or solid supports (e.g., beads) conjugated to a tag. In some cases, a partition within the set of individual partitions contains, on average, one tag or solid support (e.g., bead) conjugated to a tag. In some cases, a partition within the set of individual partitions contains, on average, 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, 30, 40, 50, 60, 70, 80, 90, or 100 tags or solid supports (e.g., beads) conjugated to a tag. In some cases, the set of individual partitions comprises empty partitions. Often, the set of individual partitions comprises some empty partitions and some partitions comprising a tag (e.g., at most one tag, at most two tags, etc.). In some cases, a tag comprises a plurality of components (e.g., a plurality of molecules). In some cases, a technique is applied to ensure that all of the partitions comprise at most one tag or bead; for example, the empty partitions (e.g., droplets) may be sorted out by a flow sorter.

Labels or tags may be added to the contents of individual partitions in order to enable later identification of the particular analyte (e.g., cell) that is the source of a particular component (e.g., molecule). In some cases, the labels or tags are partitioned into a set of partitions; the partitions may then be used to label the partitioned analytes (e.g., cells).

The labels or tags may be conjugated to beads, for example, which are partitioned into the set of partitions (e.g., droplets). In some cases, material for generating labels or tags is partitioned into a set of partitions; the tags are then produced in the individual partitions and later used to label the analytes. For example, DNA tags (e.g., free DNA tags, DNA tags conjugated to beads) may be partitioned into individual partitions and then subsequently subjected to a polymerase chain reaction (PCR) to produce copies of the tags in solution within the partitions. The tags in solutions may then be used to directly label an analyte, or as a template that is used to label the analyte in a subsequent reaction. In some cases, the DNA tags (are subjected to an in vitro transcription reaction in order to produce RNA tags within the partitions. The RNA tags in solutions may then be used to directly label an analyte, or as a template that is used to label the analyte in a subsequent reaction such as a RACE reaction.

Intra-Partition Tagging

In some cases, the contents of the analyte and tag partitions described herein can be combined in order to facilitate labeling on a per-analyte (e.g., per-cell) basis. For example, if the partitions are droplets, individual droplets from an analyte set of droplets can be merged with individual droplets from a tag set of droplets in order to facilitate labeling of the analytes. Methods of combining contents of partitions are described elsewhere herein.

In some cases, a sample containing cells can be partitioned into individual partitions, each partition comprising cell(s). Tags can be applied to the analyte components (e.g., DNA, RNA, etc.) within a partition, so that each component within a partition is labeled with the same tag, such as a cell tag capable of identifying a particular cell. In some cases, tags can be applied to the analyte components (e.g., DNA, RNA, etc.) within a partition, so that a portion of the components within a partition is labeled with the same tag, such as a cell tag capable of identifying a particular cell. The tags may comprise a molecule tag, where each tag within a partition comprises a different molecule tag. In some cases, an individual tag may comprise both a cell tag and a molecule tag. The components within a partition may be labeled with tags, so that each analyte or analyte component in a partition is labeled with a different molecule tag label. The components within a partition may be labeled with tags, so that each analyte or analyte component in a partition is labeled with an identical cell tag. The components within a partition may be labeled with tags, so that each analyte or analyte component in a partition is labeled with a different molecule tag and an identical cell tag.

FIG. 8 depicts exemplary methods provided herein. In some cases, a solution comprising beads conjugated to nucleic acids [820] is partitioned into droplets [810]. In some cases, the droplets contain at most one bead. The nucleic acids may be double-stranded DNA (dsDNA), single-stranded DNA (ssDNA), RNA, or a combination DNA and RNA. The nucleic acids may comprise a unique tag sequence (e.g. a 9-mer randomer). By “randomer” as used herein it is meant an oligonucleotide with randomly synthesized bases. When the beads are conjugated to dsDNA, in vitro transcription [801] may then be used to produce RNA transcripts from the tag dsDNA conjugated to the beads [830]. The oligonucleotides can be template switching oligonucleotides. In some cases, a solution comprising cells [850] is partitioned into droplets [840]. In some cases, the droplets contain at most one cell. Cells may be lysed [802]. Lysis can be performed by adding a lysis solution, buffer, or detergent to the cell, by sonication, by shear, by freezing and thawing, by heating, by electrical lysis, by grinding, or by any other appropriate method. Addition of lysis solution, buffer, or detergent can be performed by droplet injection, by droplet merging, or by any other appropriate method.

Droplets containing cell lysate [860] may then be merged [803] with droplets containing amplified oligonucleotides. Alternatively, the material within the droplets may be lysed following merging. Droplet merging may be conducted by any appropriate means, including passive droplet merging (e.g. at a microfluidic junction) and active droplet merging (e.g. electric, magnetic, thermal, or optical means). Alternatively, cell lysis may be conducted after droplet merging. Reverse transcription [804] with oligonucleotides and cell lysate may then be performed in the merged droplet [870] to produce cDNA [880] from a single cell tagged with the same unique tag sequence.

In some cases, oligonucleotides derived from many cells may be tagged in parallel, for example on a high-throughput basis. A solution comprising beads with surface-bound nucleic acids and a solution comprising cells may each be partitioned into droplets. In some cases, the bead droplets ([911], [921], [931]) comprise a bead ([912], [922], [932]) and the cell droplets ([914], [924], [934]) comprise a cell ([915], [925], [935]). The nucleic acids may be double-stranded DNA, single-stranded DNA, RNA, or a combination DNA and RNA. The nucleic acids may comprise a unique tag sequence (e.g. a 9-mer randomer). In vitro transcription may then be used to amplify oligonucleotides with a tag sequence ([913], [923], [933]). The oligonucleotides can be template switching oligonucleotides. Cells may be lysed. Lysis can be performed by adding a lysis solution, buffer, or detergent to the cell, by sonication, by shear, by freezing and thawing, by heating, by electrical lysis, by grinding, or by any other appropriate method. Addition of lysis solution, buffer, or detergent can be performed by droplet injection, by droplet merging, or by any other appropriate method. Droplets containing lysate ([916], [926], [936]) may then be merged with droplets containing amplified oligonucleotides. Droplet merging may be conducted by any appropriate means, including passive droplet merging (e.g. at a microfluidic junction) and active droplet merging (e.g. electric, magnetic, thermal, or optical means). Alternatively, cell lysis may be conducted after droplet merging. Reverse transcription with oligonucleotides and cell lysate may then be performed in the merged droplets ([917], [927], [937]) to produce cDNA ([918], [928], [938]), where the cDNA in each droplet is from a single cell and is tagged with the same unique tag sequence. After labeling, the cDNA from all cells may be pooled [901]. The number of cells may be at least 10, at least 100, at least 1000, at least 10,000, or at least 100,000. PCR may then be conducted on the pooled cDNA [902]. The PCR may specifically target genes or regions of interest for sequencing. The genes or regions of interest may comprise immunoglobulin heavy chain (IgH), immunoglobulin light chain (IgL), T-cell receptor beta (TCRb), T-cell receptor-alpha (TCRa), or immune cell markers. Amplified DNA may then be sequenced [903]. Information from sequencing may then be demultiplexed [904] based on the tag sequences.

The methods provided herein may result in the production of cDNA constructs which contain genetic information from cell mRNA and tag information from tags. The tag information may comprise cell tag information. The tag information may comprise molecule tag information.

Partitions may be merged into one total population subsequent to the generation of labeled molecules within the partitions. After merging, any desired analysis may be performed on the bulk population of molecules without loss of information resolution.

In some cases, a total population of cDNA molecules comprising mRNA information, cell tag information, and molecule tag information may be amplified by universal PCR to increase the total number of cDNA constructs. Sequence specific PCR may then be performed to produce sequencing-compatible DNA molecules containing cell tag information, molecule information, and genetic information of interest. By sequencing these molecules, genetic information may then be traced backwards and associated with original cells and molecules by use of tag information.

III. Tags

The methods, compositions, and kits described herein include the use of tags to identify individual sub-populations of a sample, such as cells or molecules originating from individual cells, (cell tags) and specific analytes, such as oligonucleotide sequences, contained within the individual sub-populations (molecule tags). In general, tags may comprise oligonucleotides, DNA, RNA, polypeptides, antibodies, and/or other proteins. In particular, this disclosure describes tags comprised of oligonucleotides. In some cases, the oligonucleotide comprises cell tags. In some cases, the oligonucleotides comprises molecule tags. In some cases, an oligonucleotides comprises both molecule tags and cell tags. In some cases, the oligonucleotides comprises additional sequences.

Cell Tags

Cell tags may be unique N-mer sequences that are used to identify the individual sub-population, such as the cell, from which a given component, such as gene expression product or oligonucleotide sequence (mRNA, DNA), was derived. Cell tags may comprise random N-mer oligonucleotide sequences (one unique identifier sequence per cell) that may be incorporated into cDNA constructs that include promoter sequences, adaptor sequences, and/or primer sequences that enable downstream amplification and sequencing of all or part of the gene products or oligonucleotide sequences that are being targeted. Each specific N-mer sequence may serve as a unique identifier for an individual cell, and may be incorporated into cDNA constructs using any of the several methods described herein.

The set of random N-mers used as cell tags in the examples disclosed herein are 9 bases long, but this length may be varied if a larger (or smaller) number of unique identifiers is required. In general, the length of the random N-mer used as a cell tag may range from 2 bases to 100 bases long, or may be more than 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 20, 30, 40, 50, 60, 70, 80, 90, 100, 110, 120, 130, 140, 150, 160, 170, 180, 190, or 200 bases long. Preferentially, the random N-mer used as a cell tag is at least 4 bases long and at most 14 bases long. For a random 9-mer, there are 2.62×10⁵ possible base combinations available for tagging individual cells. Examples of combinatorial properties for a given randomer length are provided in FIG. 17. In some examples presented below, the cell tag sequence initially consists of a DNA sequence, and is subsequently converted to a complementary RNA sequence incorporated into a template switching oligonucleotide using an in vitro transcription reaction. In other examples, the cell tag sequence could be introduced directly as RNA, for example, as part of a library of template switching oligonucleotides synthesized using combinatorial solid-state synthesis techniques. In many cases, the cell tag serves as both an identifier for the cell in which a given sequence is expressed, and as part of a template switching primer that is used to incorporate a collection of molecule tags and other promoter, primer, and/or adaptor sequences into the final cDNA library.

Molecule Tags

Molecule tags can comprise unique N-mer sequences that are used to identify the individual molecules (gene sequences, oligonucleotides, mRNA, DNA) that were derived from a given individual cell. Molecule tags can comprise random N-mer sequences (one unique identifier per molecule) that may be incorporated into cDNA constructs that include promoter sequences, adaptor sequences, and/or primer sequences that enable downstream amplification and sequencing of all or part of the gene products or oligonucleotide sequences that are being targeted. Each specific N-mer sequence may serve as a unique identifier for the target sequence to be detected, and is incorporated into cDNA constructs using any of the several methods described below. The set of random N-mers used as molecule tags in some of the examples disclosed herein are 9 bases long, but again, may be adjusted in length to provide the number of unique molecule identifiers required for a given application. In general, the length of the random N-mer used as a molecule tag may range from 2 bases to 100 bases long, or may be more than 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 20, 30, 40, 50, 60, 70, 80, 90, 100, 110, 120, 130, 140, 150, 160, 170, 180, 190, or 200 bases long. Preferentially, the random N-mer used as a molecule tag is at least 4 bases long and at most 14 bases long. Depending on the specific set of tagging reactions performed, the molecule tag may consist of a DNA sequence, an RNA sequence, or may be converted from one to another during the course of the tagging procedure by means of performing in vitro transcription or reverse transcription reactions. In some cases, the molecule tag serves as both an identifier for the target sequence of interest, and as part of a template switching primer that is used to incorporate other promoter, primer, and/or adaptor sequences into the final cDNA library.

Additional Oligonucleotide Sequences

In some cases, the additional nucleotide sequences comprise primers, such as PCR primers or reverse transcription primers or transcription promoter sequences. The oligonucleotide constructs used to incorporate cell and molecular tags into the cDNA library typically include primer sequences for the purpose of initiating in vitro reverse transcription reactions. Any of a number of known primer sequences may be used, for example, sequence-specific primers, oligo(dT) primers, random hexamer primers, or random decamer primers. In some cases, a general reverse transcriptase primer, for example, a poly-dT primer that complements the 3′ poly-A tail found on eukaryotic mRNA, is used to reverse transcribe all mRNA expressed in a cell. In other cases, the reverse transcriptase primer may be designed for specific gene sequences in order to perform sequence-specific reverse transcription.

In some cases, the oligonucleotides may comprise universal primer sequences. The oligonucleotide constructs used to incorporate cell and molecule tags into the cDNA library may include universal primer sequences for the purpose of initiating PCR amplification of the oligonucleotide construct. Any of a number of known universal primer sequences may be used, including the T7 universal primer sequence, or the SP6 universal primer sequence.

In some cases, the oligonucleotides may comprise promoter sequences. The oligonucleotide constructs used to incorporate cell and molecular tags into the cDNA library may also include promoter sequences for the purpose of initiating polymerase reactions. Any of a number of well-known promoter sequences, for example, the T3 promoter sequence, or the T7 promoter sequence may be used.

The oligonucleotide constructs used to incorporate cell and molecular tags into the cDNA library may also include adaptor sequences, which can be used, for example, for the purpose of facilitating sequencing of the final cDNA library on commercial sequencing platforms. Any of a number of known adaptor sequences may be used, e.g. those recommended by Illumina, Life Technologies, Pacific Biosciences, or others.

Attachment of Oligonucleotides onto a Solid Support

A variety of techniques may be used to incorporate cell and molecule tags onto a solid support.

In one approach, a single-stranded primer sequence (for example, the T7 promoter sequence) is covalently attached to a bead, and a PCR reaction using a complementary cell tag construct is performed within a compartment containing the bead to create double-stranded DNA template on the bead. The double-stranded DNA on the bead then serves as the template for subsequent in vitro transcription and PCR reactions.

In an alternative approach, a single-stranded primer sequence that incorporates a molecule tag (for example, a random N-mer sequence) is covalently attached to the bead, and a PCR reaction using a complementary cell tag construct is performed within a compartment containing the bead to create double-stranded DNA template on the bead. Again, the double-stranded DNA on the bead then serves as the template for subsequent in vitro transcription and PCR reactions.

In some approaches, the initial primer is covalently attached to the bottom surface of a microwell, and the PCR reaction using a complementary cell tag construct is performed within the well to create double-stranded DNA template molecules attached to the bottom of the well. In yet other approaches, the initial primer is covalently attached to a microarray substrate, and the PCR reaction using a complementary cell tag construct is performed within a droplet or layer of liquid in contact with the microarray to create double-stranded DNA template molecules attached to the microarray substrate. In some cases, the PCR reaction using the appropriate primers may be performed in solution, for example, within the confines of a droplet of liquid or within a microwell.

Solid Supports

Suitable solid phase carriers include, but are not limited to, other particles, fibers, beads and or supports which have an affinity for DNA, RNA, double-stranded DNA, single stranded DNA, ssRNA and which can embody a variety of shapes, that are either regular or irregular in form, provided that the shape maximizes the surface area of the solid phase, and embodies a carrier which is amenable to microscale manipulations

An example of a suitable bead is a porous or nonporous polymer bead comprised of a copolymer of vinyl aromatic monomers. Examples of vinyl aromatic monomers include styrene, alkyl substituted styrene, alpha-methylstyrene, and alkyl substituted alpha-methylstyrene.

Another example of a suitable bead is a porous or nonporous particle such as silica, silica carbide, silica nitrite, titanium oxide, aluminum oxide, zirconium oxide which is modified to have a hydrophobic surface.

In some aspects the solid surface has a functionalize surface such that the surface is coated with moieties which reversibly bind nucleic acid (e.g., DNA, RNA). One example is a surface which is coated with moieties which each have a free functional group which is bound to the amino group of the amino silane or the microparticle; as a result, the surfaces of the microparticles are coated with the functional group containing moieties. The functional group acts as a bioaffinity adsorbent for polyalkylene glycol precipitated DNA. In one embodiment, the functional group is a carboxylic acid. A suitable moiety with a free carboxylic acid functional group is a succinic acid moiety in which one of the carboxylic acid groups is bonded to the amine of amino silanes through an amide bond and the second carboxylic acid is unbonded, resulting in a free carboxylic acid group attached or tethered to the surface. Suitable solid phase carriers having a functional group coated surface that reversibly binds nucleic acid molecules are for example, magnetically responsive solid phase carriers having a functional group-coated surface, such as, but not limited to, amino-coated, carboxyl-coated and encapsulated carboxyl group-coated paramagnetic microparticles.

IV. Template Switching Oligonucleotides

Tags may be delivered by template switching oligonucleotides. Template switching oligonucleotides may be RNA or DNA. Template switching oligonucleotides may comprise a hybridization region and a template region. The hybridization region can comprise any sequence capable of hybridizing to the target. In some cases, the hybridization region comprises a series of G bases to complement the overhanging C bases at the 3′ end of a cDNA molecule, which result from the terminal transferase activity of reverse transcriptase. The series of G bases may comprise 1 G base, 2 G bases, 3 G bases, 4 G bases, 5 G bases, or more than 5 G bases. The template region can comprise any sequence to be incorporated into the cDNA by primer extension reaction. In some cases, the template region comprises one or more tag sequences. The template region can contain one tag sequence, two tag sequences, three tag sequences, or more than three tag sequences. The template region may comprise a molecule tag region and a cell tag region.

Template switching oligonucleotides may further comprise additional regions. In some cases, the template switching oligonucleotide comprises one or more sequencing adaptor sequences or partial sequencing adaptor sequences. The sequencing adaptor sequence may be an Illumina adaptor sequence. The sequencing adaptor sequence may be a 454 adaptor sequence. The sequencing adaptor sequence may be a SOLiD adaptor sequence. In some cases, the template switching oligonucleotide comprises one or more universal sequences. Universal sequences may be used for primer binding for universal PCR or other reactions.

The template switching oligonucleotide (e.g., RNA oligonucleotide) may be incorporated into the cDNA sequence during the extension reaction of the reverse transcription. The hybridization region allows the template switching oligonucleotide to hybridize to the cDNA. Once the reverse transcriptase enzyme reaches the end of the mRNA, it switches to the template switching oligonucleotide and incorporates its sequence into the cDNA as well. See, e.g., FIG. 3.

The assays provided herein can comprise a Rapid Amplification of cDNA Ends (RACE) assay as depicted, for example, in FIGS. 3-5. In some cases, the RACE assay may be performed on mRNA transcripts [310] from a cell. The mRNA transcript can be from a variety of genes, including genes related to immune function, such as IgH or TCR beta. In some cases, the mRNA transcript comprises six regions: a 5′ untranslated region (5′ UTR), a leader region, a V region, a D region, a J region, and a C region. Reverse transcription (RT) primers [320] may then be added [301]. In some cases, the RT primers can be specific for the heavy chain constant region of an immunoglobulin, e.g. IgG, IgA, IgM, IgE, or IgD, kappa (IgK) and lambda (IgL)light chain regions of a immunoglobulin, T cell receptor alpha constant regions (TRAC) and or T cell receptor beta constant regions(TRBC). Reverse transcriptase can then be used to perform reverse transcription [302] of the mRNA, resulting in cDNA [330] with overhanging C bases from the terminal transferase activity of the reverse transcriptase.

Template switching oligonucleotides [340] may then be added, switching the template molecule from the original mRNA to the cDNA. In some cases, the template switching oligonucleotide comprises three regions: a partial sequencing adaptor sequence (i.e., a universal sequence to the 3′ end), a randomer molecule tag sequence, and a homopolymerice tail (i.e., G base region to complement the overhanging C bases of the cDNA). The partial sequencing adaptor sequence may be a partial Illumina adaptor sequence. The randomer molecule tag sequence may consist of nine random bases (9mer). The cDNA can then be extended [303] to include the template switching oligonucleotide sequence, resulting in a full length cDNA [350]. In some cases, the full length cDNA comprises a unique randomer molecule tag sequence at the V end and a sequence complementary to part of a sequencing adaptor sequence. The sequencing adaptor sequence can be an Illumina adaptor sequence.

Additional primers may then be added [401]. In some cases, three additional primers are added, comprising: nested C primers with a partial sequencing adaptor sequence [420], a sequencing adaptor sequence with a 6 base pair index [430] referred to as a long primer, and a short primer [440]. The sequencing adaptor sequences can be Illumina adaptor sequences. PCR may then be conducted [402] to produce amplified cDNA [450] which incorporates the primer sequences. In some cases, the PCR is conducted for 30 cycles. Additional primers may then be added. In some cases, one additional primer is added, comprising a full-length sequencing adaptor sequence [460]. The sequencing adaptor sequence can be an Illumina adaptor sequence, such as P7. PCR may then be conducted [403] to add on the full-length adaptor sequence. In some cases, the PCR is conducted for 8 cycles.

The resulting cDNA products [510] may then be sequenced and analyzed. In some cases, this is accomplished by sequencing. The sequencing may be performed with the use of read primers. In some cases, three read primers are used, comprising: a first read primer (Read1) [520], with a 150 base pair read area [530] covering most of the CDR3 region, a second read primer (Read2) [540], with a 150 base pair read area [550] covering the unique randomer molecule tag barcode sequence, the 5′UTR, the leader region, and part of the V region, and an index read primer [560], with a read area covering the index region. In some aspects the read length is 50, 100, 150, 200, 250, 500, 750, 1000, 1250 or 1500 base pairs. The cDNA products may be analyzed by other techniques, including hybridization probes.

The RACE assay may be 5′ RACE or 3′ RACE. In some cases, the sequencing is performed from linear cDNA products. In some cases, the sequencing is performed from circularized cDNA products. Circularized cDNA products may be produced by melting and circularizing, that is, the 3′ and 5′ ends of the linear DNA are brought together and bonded (e.g. by a circular ligase enzyme such as CircLigase™, EPICENTRE® Biotechnologies).

V. Reactions

The methods and compositions provided in this disclosure can comprise one or more reactions involved in performing analysis. The reactions can comprise sample preparation such as cell lysis, tag production such as in vitro transcription, tagging such as reverse transcription with template switching oligonucleotides, analyte amplification such as PCR, and analysis such as sequencing.

Cell Lysis

Following partitioning of cells into individual partitions, the cells may be lysed to release the intracellular contents, including RNA, DNA, proteins, and other intracellular components. The lysis reaction may be performed using any number of well-known techniques, for example by the addition of a detergent, addition of a lysis solution or hypotonic buffer, use of mechanical agitation or sonication, repeated cycles of freezing and thawing, or by any other appropriate method. Addition of lysis solutions, buffers, or detergents can be performed by conventional liquid dispensing methods, by droplet injection or merging techniques, or by any other appropriate method. Partitions containing cell lysate [860] may subsequently be merged [803] with partitions containing tags prior to performing downstream reactions. In some cases, lysis can be performed prior to merging the contents of the partition with another partition, such as a partition containing tag. In some cases, lysis can be performed after merging the contents of a partition with other partitions, such as partitions containing tags.

In Vitro Transcription.

In vitro transcription reactions may be used to convert DNA template molecules into RNA transcripts. The DNA template molecules used in the methods and compositions disclosed herein may be designed to contain promoter sequences for facilitating in vitro transcription. The DNA template molecules may comprise random or unique N-mer tagging sequences, and may also include primer and adaptor sequences for facilitating downstream amplification and sequencing reactions. The products of in vitro transcription reactions in the present method are RNA template switching oligonucleotides comprised of the said appropriate primer, adaptor, and/or tag sequences. In some cases, double-stranded DNA templates bound to beads are transcribed to produce RNA constructs comprising template switching oligonucleotides, which may include primer, adaptor, and/or random or unique N-mer cell tag sequences, random or unique N-mer molecule tag sequences, or both random or unique N-mer cell tag and N-mer molecule tag sequences. Alternatively, in vitro transcription may be used with free template DNA to create the RNA template switching constructs in solution.

In some cases, in vitro transcription reactions may be performed prior to combining DNA template sequences with cell lysates or other biological samples. In some cases, in vitro transcription reactions may be performed after combining DNA template sequences with cells, or material derived from a cell (e.g., polynucleotides, polypeptides, etc.), or other biological samples. DNA templates used in the in vitro transcription reaction can comprise a promoter sequence for facilitating the transcription reaction. Any number of known promoter sequences may be used, for example, the T3 RNA polymerase promoter sequence, T7 RNA polymerase promoter sequence, or any other suitable RNA polymerase promoter sequence. The transcription reaction may be performed under the appropriate reaction conditions using any appropriate in vitro transcription system, for example, the T3 RNA polymerase system, which includes the appropriate transcriptase (T3 RNA polymerase), ribonucleotide triphosphates, and buffer components (for example, dithiothreitol and magnesium ions).

The in vitro transcription assays can be designed so that all of the oligonucleotides present in a sample are transcribed into RNA, or only a portion of the oligonucleotides are transcribed into RNA. For example, all of the oligonucleotides may contain the same promoter; or only a subset of the oligonucleotides may contain the promoter. In some cases, the oligonucleotides comprise different promoters.

Polymerase Chain Reaction (PCR)

PCR reactions may be used at any step described herein. PCR reactions may be used, for example, to (i) perform bead-based synthesis of DNA constructs that incorporate promoters, primers, adaptors, and/or unique N-mer cell tags, (ii) perform bead-based synthesis of DNA constructs that incorporate promoters, primers, adaptors, and/or unique N-mer molecule tags, (iii) perform bead-based synthesis of DNA constructs that incorporate promoters, primers, adaptors, and/or both unique N-mer cell tags and/or unique N-mer molecule tags, (iv) perform general amplification of pooled cDNA libraries, or (v) to perform sequence-specific amplification of cDNA sequences that represent the target genes or oligonucleotide sequences of interest. In some preferred approaches, emulsion PCR using either the AmpliTaq or Phusion polymerase, single-stranded primer (e.g. Roche 454 primer A) attached to beads (e.g. Roche 454 beads), and a DNA construct comprising the 454 A and/or 454 B primer sequences, The T7 promoter, a partial Illumina adaptor sequence, an N-mer cell tag, a universal adapter sequence ending in GGG, and an optional restriction site sequence is used to synthesize double-stranded DNA on beads, which are then used to perform bead-based in vitro transcription (following treatment of the beads with a blunt-end restriction digest, if necessary to remove the restriction site) to produce RNA template switching oligos that incorporate the N-mer cell tag.

Although many of the examples disclosed herein describe bead-based synthesis of DNA constructs using emulsion PCR, the PCR synthesis and amplification reactions may also be performed in solution to synthesize said DNA constructs. Any of a number of known DNA polymerase systems (including the polymerase, suitable primers, necessary cofactors and buffers, and deoxyribonucleotide triphosphates) may be used under the appropriate set of reaction conditions (e.g., suitable melt temperatures, annealing temperatures, and elongation reaction temperatures) to perform the synthesis or amplification. Exemplary systems include the Phusion polymerase system, the Taq polymerase system, the AmpliTaq polymerase system, or any other suitable polymerase system. Preferentially, the polymerase chosen for use in PCR reactions used to synthesize the DNA constructs described above is one that does not add overhanging A bases at the 3′ end of the construct, as this necessitates removal using T4 DNA polymerase or other suitable exonuclease.

In some cases, PCR reactions are used with an appropriate set of primers in a general fashion to synthesize and/or amplify all DNA template or cDNA library molecules. In other cases, PCR reactions are used with an appropriate set of primers in a selective fashion to specifically amplify those target gene sequences or target oligonucleotide sequences present in the cDNA library. In some cases, both general and specific PCR reactions are performed. The number of PCR reaction cycles utilized may vary from about 2 cycles to 40 cycles or more.

DNA Sequencing

DNA sequencing may be used to sequence the pooled and amplified cDNA library in order to identify: (i) the identity of the cell tag, (ii) the identity of the molecule tag, and/or (iii) the complete or partial sequence of the target genes or oligonucleotides of interest. DNA sequencing is performed using any of a number of commercially-available sequencing systems (i.e. reagents, kits, and instruments), for example the lllumina MiSeq, HiSeq, or NextSeq 500 systems; the Life Technologies SOLID sequencing system, the Pacific Biosciences SMRT sequencing system, or any of the other commercially-available or emerging sequencing technology platforms. The choice of adapter sequences incorporated into the design of the DNA or RNA constructs used for adding cell tags and/or molecule tags to the cDNA library is driven by the choice of sequencing system used.

In certain embodiments, the sequencing technique used in the methods of the provided invention generates at least 100 reads per run, at least 200 reads per run, at least 300 reads per run, at least 400 reads per run, at least 500 reads per run, at least 600 reads per run, at least 700 reads per run, at least 800 reads per run, at least 900 reads per run, at least 1000 reads per run, at least 5,000 reads per run, at least 10,000 reads per run, at least 50,000 reads per run, at least 100,000 reads per run, at least 500,000 reads per run, at least 1,000,000 reads per run, at least 2,000,000 reads per run, at least 3,000,000 reads per run, at least 4,000,000 reads per run at least 5000,000 reads per run, at least 6,000,000 reads per run at least 7,000,000 reads per run at least 8,000,000 reads per run, at least 9,000,000 reads per run, or at least 10,000,000 reads per run.

In some embodiments the number of sequencing reads per B cell sampled should be at least 2 times the number of B cells sampled, at least 3 times the number of B cells sampled, at least 5 times the number of B cells sampled, at least 6 times the number of B cells sampled, at least 7 times the number of B cells sampled, at least 8 times the number of B cells sampled, at least 9 times the number of B cells sampled, or at least at least 10 times the number of B cells The read depth allows for accurate coverage of B cells sampled, facilitates error correction, and ensures that the sequencing of the library has been saturated.

In some embodiments the number of sequencing reads per T-cell sampled should be at least 2 times the number of T-cells sampled, at least 3 times the number of T-cells sampled, at least 5 times the number of T-cells sampled, at least 6 times the number of T-cells sampled, at least 7 times the number of T-cells sampled, at least 8 times the number of T-cells sampled, at least 9 times the number of T-cells sampled, or at least at least 10 times the number of T-cells The read depth allows for accurate coverage of T-cells sampled, facilitates error correction, and ensures that the sequencing of the library has been saturated.

VI. Mechanisms for Partitioning and Merging

The methods and compositions provided in this disclosure may be implemented in partitions. These partitions may contain single cells for analysis of a sample at the single cell level. These partitions may contain tags or means for delivering tags. These partitions may comprise droplets, microwell arrays, microarrays, or any other suitable technology.

A. Droplets

In some cases, partitioning is performed with droplets. Droplets may comprise an aqueous medium (e.g. water, buffer solution, cell growth medium) surrounded by an immiscible oil phase (e.g. mineral oil, silicone oil, perfluorinated oil). Droplets may comprise an aqueous medium (e.g. water, buffer solution, cell growth medium) surrounded by air. Droplets may be located in a vial, tube, capillary, syringe, microfluidic channel, or on a surface (e.g. a piezoelectric surface for droplet actuation by surface acoustic waves). Droplets may be generated by any appropriate method, such as a microfluidic device (e.g., RainDance RainDrop system, microfluidic T-junction), bulk emulsification, or pipetting.

Droplets may additionally be formed with surfactants. Examples of surfactants include Triton X-100, SDS, ABIL EM90, Span80, monolein, oleic acid, Tween 20, Tween 80, Synperonic PEF, C12E8, n-Butanol, phospholipids, PF-octanol, PF-decanol, PF-TD OEG, PFPE-COOH, PFPE-COONH4, PFPE-PEG, PFPE-DMP, and Pico-Surf.

Droplets may be sorted. For example, a sample containing cells may be partitioned into droplets containing cells (e.g., at most one cell) and droplets that contain no cells. Droplets containing zero cells may be sorted and discarded. Droplets may be sorted by any appropriate method, such as a microfluidic droplet sorting device or fluorescence-activated droplet sorting (FADS) system.

Droplets may be merged. For example, a droplet containing a cell or cell lysate may be merged with a droplet containing tags. Droplet merging may be conducted by any appropriate means, including passive droplet merging (e.g. at a microfluidic junction) and active droplet merging (e.g. electric, magnetic, thermal, or optical means).

The number of cells partitioned may be at least 10, at least 100, at least 1000, at least 10,000, or at least 100,000. Partitions may include at most one cell, at most two cells, at most three cells, at most four cells, at most five cells, or at most ten cells.

B. Microwell Arrays

In some cases, partitioning is performed with one or more microwell arrays.

Microwell arrays may comprise plates, films, tapes, or other substrates comprising an array of microwells. The substrate may comprise plastic (e.g. polystyrene, polypropylene, polycarbonate, cyclo-olefin, Lucite), epoxy, photoresist (e.g. SU-8), PDMS, glass, metal, or any other suitable material. Microwell arrays may comprise at least 1, at least 6, at least 24, at least 96, at least 384, at least 1536, at least 3456, or at least 9600 wells. Sample may be distributed in the microwells by pipetting, liquid handlers, robotic fluid handlers, or any other appropriate method.

Microwell contents may be sorted. For example, a sample containing cells may be partitioned into microwells containing at most one cell, according to Poisson distribution statistics. Microwell contents containing zero cells may be chosen to be discarded. Microwell contents may be sorted by any appropriate method, such as fluorescence or optical detection. Microwell contents to be discarded may be removed from the microwells by pipetting, liquid handlers, robotic fluid handlers, or any other appropriate method.

Microwell contents may be merged. For example, the contents of a microwell with a cell or cell lysate may be merged with the contents of a microwell containing tags. Microwell content merging may be conducted by any appropriate method, including pipetting, liquid handling, robotic fluid handling, or physical alignment and contacting of multiple microwell plates.

The number of cells partitioned may be at least 10, at least 100, at least 1000, at least 10,000, or at least 100,000. Partitions may include at most one cell, at most two cells, at most three cells, at most four cells, at most five cells, or at most ten cells. In some cases, the microwells are designed to accommodate not more than one cell and/or not more than one bead. For example, the diameter of the microwells may be less than 10, 15, 20, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, 100, 105, 100, 110, 120, 130, 140, 150, 160, 170, 180, 190, 200, 250, or 300 um.

C. Microarrays with Dividers

In some cases, partitioning is performed with one or more microarray chips. Microarray chips may comprise chips, slides, or other substrates comprising an array of locations with bound material. The substrate may comprise glass, silicon, nylon, plastic, or any other suitable material. The material bound to the array locations may comprise DNA, RNA, proteins, peptides, antibodies, cells, chemicals, tissues, or carbohydrates.

Material bound to microarray locations may be further partitioned by inclusion of dividers or walls. A membrane or material layer comprising dividers or walls may be aligned and bonded to the microarray surface, partitioning each microarray location and providing a partitioned volume. The membrane may comprise PDMS, glass, silicon, plastic (e.g. PMMA), metal, or any other suitable material.

Subsequent to partitioning by divider, the contents of a microarray location may be detached from the substrate, or may be replicated, to produce unbound copies of the location contents. For example, a microarray location containing surface-bound DNA may be reacted to produce RNA molecules in free solution by in vitro transcription.

Microarray location volume contents may be merged. For example, the contents of a microarray location volume containing tags may be merged with a volume containing cells or cell lysate. Microarray location volume content merging may be conducted by any appropriate method, including pipetting, liquid handling, robotic fluid handling, or physical alignment and contacting with microwell plates or other partitioned microarrays.

VI. Digital Counting and Analysis

The high throughput sequencing provides a very large dataset.

High-throughput analysis can be achieved using one or more bioinformatics tools, such as ALLPATHS (a whole genome shotgun assembler that can generate high quality assemblies from short reads), Arachne (a tool for assembling genome sequences from whole genome shotgun reads, mostly in forward and reverse pairs obtained by sequencing cloned ends, BACCardl (a graphical tool for the validation of genomic assemblies, assisting genome finishing and intergenome comparison), CCRaVAT & QuTie (enables analysis of rare variants in large-scale case control and quantitative trait association studies), CNV-seq (a method to detect copy number variation using high throughput sequencing), Elvira (a set of tools/procedures for high throughput assembly of small genomes (e.g., viruses)), Glimmer (a system for finding genes in microbial DNA, especially the genomes of bacteria, archaea and viruses), gnumap (a program designed to accurately map sequence data obtained from next-generation sequencing machines), Goseq (an R library for performing Gene Ontology and other category based tests on RNA-seq data which corrects for selection bias), ICAtools (a set of programs useful for medium to large scale sequencing projects), LOCAS, a program for assembling short reads of second generation sequencing technology, Maq (builds assembly by mapping short reads to reference sequences, MEME (motif-based sequence analysis tools, NGSView (allows for visualization and manipulation of millions of sequences simultaneously on a desktop computer, through a graphical interface, OSLay (Optimal Syntenic Layout of Unfinished Assemblies), Perm (efficient mapping for short sequencing reads with periodic full sensitive spaced seeds, Projector (automatic contig mapping for gap closure purposes), Qpalma (an alignment tool targeted to align spliced reads produced by sequencing platforms such as Illumina, Solexa, or 454), RazerS (fast read mapping with sensitivity control), SHARCGS (SHort read Assembler based on Robust Contig extension for Genome Sequencing; a DNA assembly program designed for de novo assembly of 25-40 mer input fragments and deep sequence coverage), Tablet (next generation sequence assembly visualization), and Velvet (sequence assembler for very short reads).

A Non-limiting example of data analysis steps are summarized below:

Grouping reads with the same cellular and/or molecular tag: Initially sequences are matched based on identical cellular and or molecular tags.

Build a minimum spanning forest for each group: Cluster into sungroups (trees) if Hamming distance is greater than 5%.

For each subgroup (or tree), create a vector of sums of correct probabilities for each called base in each read.

Construct a consensus read from the base with the maximum sum in each position: Consensus reads are used for mutation analysis and diversity measurement.

VII. Information Systems

The methods provided herein may be processed on a server or a computer server (FIG. 18). The server [1801] includes a central processing unit (CPU, also “processor”) [1805] which can be a single core processor, a multi core processor, or plurality of processors for parallel processing. A processor used as part of a control assembly may be a microprocessor. The server [1801] also includes memory [1810] (e.g. random access memory, read-only memory, flash memory); electronic storage unit [1815] (e.g. hard disk); communications interface [1820] (e.g. network adaptor) for communicating with one or more other systems; and peripheral devices [1825] which may include cache, other memory, data storage, and/or electronic display adaptors. The memory [1810], storage unit [1815], interface [1820], and peripheral devices [1825] are in communication with the processor [1805] through a communications bus (solid lines), such as a motherboard. The storage unit [1815] can be a data storage unit for storing data. The server [1801] is operatively coupled to a computer network (“network”) [1830] with the aid of the communications interface [1820]. A processor with the aid of additional hardware may also be operatively coupled to a network. The network [1830] can be the Internet, an intranet and/or an extranet, an intranet and/or extranet that is in communication with the Internet, a telecommunication or data network. The network [1830] in some cases, with the aid of the server [1801], can implement a peer-to-peer network, which may enable devices coupled to the server [1801] to behave as a client or a server. In general, the server may be capable of transmitting and receiving computer-readable instructions (e.g., device/system operation protocols or parameters) or data (e.g., sensor measurements, raw data obtained from detecting nucleic acids, analysis of raw data obtained from detecting nucleic acids, interpretation of raw data obtained from detecting nucleic acids, etc.) via electronic signals transported through the network [1830]. Moreover, a network may be used, for example, to transmit or receive data across an international border.

VIII. Applications

The invention finds use in the prevention, treatment, detection, diagnosis, prognosis, or research into any condition or symptom of any condition, including cancer, inflammatory diseases, autoimmune diseases, allergies and infections of an organism (i.e., bacterial, viral or fungal). The organism is preferably a human subject but can also be derived from non-human subjects, e.g., non-human mammals. Examples of non-human mammals include, but are not limited to, non-human primates (e.g., apes, monkeys, gorillas), rodents (e.g., mice, rats), cows, pigs, sheep, horses, dogs, cats, or rabbits.

Examples of cancer include prostrate, pancreas, colon, brain, lung, breast, bone, and skin cancers. Examples of inflammatory conditions include irritable bowel syndrome, ulcerative colitis, appendicitis, tonsilitis, dermatitis. Examples of atopic conditions include allergy, asthma, etc. Examples of autoimmune diseases include IDDM, RA, MS, SLE, Crohn's disease, Graves' disease, etc. Autoimmune diseases also include Celiac disease, and dermatitis herpetiformis. For example, determination of an immune response to cancer antigens, autoantigens, pathogenic antigens, vaccine antigens, and the like is of interest.

One specific application of the claimed methods is to assess the diversity of the immune repertoire of a subject.

VDJ lineage diversity: VDJ usage is enumerated by the number of observed lineages falling into each VJ, VDJ, VJC, or VDJC (e.g., VDJ) combination at a given read-depth.

VDJ and unique sequence abundance histograms: Histograms are plotted by binning VDJ and unique sequence abundances (the latter which is either clustered or has undergone lineage-analysis filtering and grouping) into log-spaced bins.

3D representation of VJ, VDJ, VJC, or VDJC (e.g., VDJ) usage: Repertoires are represented by applying V-, D-, J-, and/or C-segments to different axes on a three-dimensional plot. Using either abundance (generally read number, which can be bias-normalized) or observed lineage diversity, bubbles of varying sizes are used at each V/D/J/C coordinate to represent the total usage of that combination.

Mutation vs. sequence abundance plots: After undergoing lineage analysis, unique sequences are binned by read-number (or bias-normalized abundance) into log-spaced bins. For a given abundance-bin, the number of mutations per unique sequence is averaged, giving a mutation vs. abundance curve.

Correlative measures of V, D, J, C, VJ, VDJ, VJC, VDJC, antibody heavy chain, antibody light chain, CDR3, or T-cell receptor usage (Pearson, KL divergence): VJ, VDJ, VJC, or VDJC (e.g., VDJ) combinations are treated as vectors with indexed components v, weighted by either lineage-diversity or abundance for that VDJ combination. Pearson correlations and KL-divergences between each pair of individuals are then calculated over the indices.

The results of the analysis may be referred to herein as an immune repertoire analysis result, which may be represented as a dataset that includes sequence information, representation of V, D, J, C, VJ, VDJ, VJC, VDJC, antibody heavy chain, antibody light chain, CDR3, or T-cell receptor usage, representation for abundance of V, D, J, C, VJ, VDJ, VJC, VDJC, antibody heavy chain, antibody light chain, CDR3, or T-cell receptor and unique sequences; representation of mutation frequency, correlative measures of VJ V, D, J, C, VJ, VDJ, VJC, VDJC, antibody heavy chain, antibody light chain, CDR3, or T-cell receptor usage, etc. Such results may then be output or stored, e.g. in a database of repertoire analyses, and may be used in comparisons with test results, reference results, and the like.

After obtaining an immune repertoire analysis result from the sample being assayed, the repertoire can be compared with a reference or control repertoire to make a diagnosis, prognosis, analysis of drug effectiveness, or other desired analysis. A reference or control repertoire may be obtained by the methods of the invention, and will be selected to be relevant for the sample of interest. A test repertoire result can be compared to a single reference/control repertoire result to obtain information regarding the immune capability and/or history of the individual from which the sample was obtained. Alternately, the obtained repertoire result can be compared to two or more different reference/control repertoire results to obtain more in-depth information regarding the characteristics of the test sample. For example, the obtained repertoire result may be compared to a positive and negative reference repertoire result to obtain confirmed information regarding whether the phenotype of interest. In another example, two “test” repertoires can also be compared with each other. In some cases, a test repertoire is compared to a reference sample and the result is then compared with a result derived from a comparison between a second test repertoire and the same reference sample.

Determination or analysis of the difference values, i.e., the difference between two repertoires can be performed using any conventional methodology, where a variety of methodologies are known to those of skill in the array art, e.g., by comparing digital images of the repertoire output, by comparing databases of usage data, etc.

A statistical analysis step can then be performed to obtain the weighted contribution of the sequence prevalence, e.g. V, D, J, C, VJ, VDJ, VJC, VDJC, antibody heavy chain, antibody light chain, CDR3, or T-cell receptor usage, mutation analysis, etc. For example, nearest shrunken centroids analysis may be applied as described in Tibshirani et at. (2002) P.N.A.S. 99:6567-6572 to compute the centroid for each class, then compute the average squared distance between a given repertoire and each centroid, normalized by the within-class standard deviation.

A statistical analysis may comprise use of a statistical metric (e.g., an entropy metric, an ecology metric, a variation of abundance metric, a species richness metric, or a species heterogeneity metric.) in order to characterize diversity of a set of immunological receptors. Methods used to characterize ecological species diversity can also be used in the present invention. See, e.g., Peet, Annu Rev. Ecol. Syst. 5:285 (1974). A statistical metric may also be used to characterize variation of abundance or heterogeneity. An example of an approach to characterize heterogeneity is based on information theory, specifically the Shannon-Weaver entropy, which summarizes the frequency distribution in a single number. See, e.g., Peet, Annu Rev. Ecol. Syst. 5:285 (1974).

The classification can be probabilistically defined, where the cut-off may be empirically derived. In one embodiment of the invention, a probability of about 0.4 can be used to distinguish between individuals exposed and not-exposed to an antigen of interest, more usually a probability of about 0.5, and can utilize a probability of about 0.6 or higher. A “high” probability can be at least about 0.75, at least about 0.7, at least about 0.6, or at least about 0.5. A “low” probability may be not more than about 0.25, not more than 0.3, or not more than 0.4. In many embodiments, the above-obtained information is employed to predict whether a host, subject or patient should be treated with a therapy of interest and to optimize the dose therein.

In addition the methods described herein are of interest as a means of detecting the earliest changes along a disease pathway (e.g., a carcinogenesis pathway, inflammatory pathway, etc.), and/or to monitor the efficacy of various therapies and preventive interventions.

The methods disclosed herein can also be utilized to analyze the effects of agents on cells. For example, analysis of changes in gene expression following exposure to one or more test compounds can performed to analyze the effect(s) of the test compounds on an individual.

Agents to be analyzed for potential therapeutic value can be any compound, small molecule, protein, lipid, carbohydrate, nucleic acid or other agent appropriate for therapeutic use. Preferably tests are performed in vivo, e.g. using an animal model, to determine effects on the immune repertoire.

Agents of interest for screening include known and unknown compounds that encompass numerous chemical classes, primarily organic molecules, which may include organometallic molecules, genetic sequences, etc. An important aspect of the invention is to evaluate candidate drugs, including toxicity testing; and the like.

Test compounds include all of the classes of molecules described above, and can further comprise samples of unknown content. Of interest are complex mixtures of naturally occurring compounds derived from natural sources such as plants, fungi, bacteria, protists or animals. While many samples will comprise compounds in solution, solid samples that can be dissolved in a suitable solvent may also be assayed. Samples of interest include environmental samples, e.g., ground water, sea water, mining waste, etc., biological samples, e.g. lysates prepared from crops, tissue samples, etc.; manufacturing samples, e.g. time course during preparation of pharmaceuticals; as well as libraries of compounds prepared for analysis; and the like (e.g., compounds being assessed for potential therapeutic value, i.e., drug candidates).

Compounds, including candidate agents, are obtained from a wide variety of sources including libraries of synthetic or natural compounds. For example, numerous means are available for random and directed synthesis of a wide variety of organic compounds, including biomolecules, including expression of randomized oligonucleotides and oligopeptides. Alternatively, libraries of natural compounds in the form of bacterial, fungal, plant and animal extracts are available or readily produced. Additionally, natural or synthetically produced libraries and compounds are readily modified through conventional chemical, physical and biochemical means, and may be used to produce combinatorial libraries. Known pharmacological agents may be subjected to directed or random chemical modifications, such as acylation, alkylation, esterification, amidification, etc. to produce structural analogs.

While preferred embodiments of the present invention have been shown and described herein, it will be obvious to those skilled in the art that such embodiments are provided by way of example only. Numerous variations, changes, and substitutions will now occur to those skilled in the art without departing from the invention. It should be understood that various alternatives to the embodiments of the invention described herein may be employed in practicing the invention. It is intended that the following claims define the scope of the invention and that methods and structures within the scope of these claims and their equivalents be covered thereby.

EXAMPLES

Example 1—Preparing Beads with Cell Tags and Molecular Tags

Beads conjugated to dsDNA comprising cellular tags may be prepared as shown in FIG. 10. A bead with DNA bound to the surface is placed in solution with oligonucleotides, primers, enzymes, and other necessary reagents for performing PCR. In some cases the DNA [1011] bound to the bead [1010] comprises one region: a promoter region. The promoter region can be a T7 promoter. In some cases, the oligonucleotides [1020] in solution comprise four regions: a promoter region [1021], a partial sequencing adaptor [1022], a randomer cell tag [1023], and a universal sequence [1024]. The promoter region can be a T7 promoter. The partial sequencing adaptor can be a partial Illumina adaptor. The universal sequence can end with GGG. In some cases, the primers in the solution [1034] comprise a universal sequence. PCR is then conducted. In some cases, the PCR is emulsion PCR [1001]. The PCR reaction may result in additional sequences being added to the DNA bound to the bead. In some cases, the resulting DNA is double-stranded DNA. In some cases the dsDNA [1040] comprises six regions: a promoter region [1041], a partial sequencing adaptor region [1042], a randomer cell tag [1043], and a universal sequence [1044]. The promoter region can be a T7 promoter. The partial sequencing adaptor region can be a partial Illumina adapter. The universal sequence can end with GGG. In vitro transcription is then conducted [1002], producing RNA oligonucleotide products [1050]. In some cases, the oligonucleotide products comprise three regions: a partial sequencing adaptor region [1052], a randomer cell tag [1053], and a universal sequence [1054]. The partial sequencing adaptor region can be a partial Illumina adapter. The universal sequence can end with GGG.

Beads conjugated to dsDNA comprising cell tags and molecule tags may be prepared as depicted in FIG. 11. A bead with DNA bound to the surface is placed in solution with oligonucleotides, primers, enzymes, and other necessary reagents for performing PCR. In some cases the DNA [1120] bound to the bead [1110] comprises four regions: a promoter region [1121], a partial sequencing adaptor region [1122], a randomer molecule tag [1123] where each DNA molecule on the bead contains a different molecule tag sequence, and a first universal sequence [1124]. The promoter region can be a T7 promoter. The partial sequencing adaptor region can be a partial Illumina adapter. In some cases, the oligonucleotides [1130] in solution comprise three regions: a first universal sequence [1134], a randomer cell tag [1135], and a second universal sequence [1136]. The second universal sequence can end with GGG. In some cases, the primers in the solution [1146] comprise a second universal sequence. PCR is then conducted. In some cases, the PCR is emulsion PCR [1101]. As a result the DNA on the bead is extended, thereby incorporating additional sequences. The reaction can also be designed so that the resulting DNA conjugated to the bead is dsDNA. In some cases the dsDNA [1150] comprises six regions: a promoter region [1151], a partial sequencing adaptor region [1152], a randomer molecule tag [1153] where each dsDNA molecule on the bead contains a different molecule tag sequence, a first universal sequence [1154], a randomer cell tag [1155], and a second universal sequence [1156]. The promoter region can be a T7 promoter. The partial sequencing adaptor region can be a partial Illumina adapter. The second universal sequence can end with GGG. In vitro transcription is then conducted [1102], producing RNA oligonucleotide products [1160]. In some cases, the oligonucleotide products comprise five regions: a partial sequencing adaptor region [1162], a randomer molecule tag [1163] where oligonucleotide products from different dsDNA molecules on the bead contain different molecule tag sequences, a first universal sequence [1164], a randomer cell tag [1165], and a second universal sequence [1166]. The partial sequencing adaptor region can be a partial Illumina adapter. The second universal sequence can end with GGG.

FIG. 12 depicts another methods for preparing beads conjugated to dsDNA comprising a cell tag. In some cases the DNA [1211] bound to the bead [1210] comprises one region: an adaptor region [1211]. The adaptor region can be a 454 A adaptor. In some cases, the oligonucleotides [1220] in solution comprise seven regions: a first adaptor region [1221], a promoter region [1222], a partial sequencing adaptor region [1223], a randomer cell tag region [1224], a universal sequence [1225], a restriction site [1226], and a second adaptor region [1227]. The first and second adaptor regions can be 454 A and 454 B adaptors, respectively. The promoter region can be a T7 promoter. The partial sequencing adaptor region can be a partial Illumina adaptor. The universal sequence can end with GGG. In some cases, the primers in the solution [1237] comprise a second adaptor region sequence. The second adaptor region can be a 454 B adaptor. PCR is then conducted. In some cases, the PCR is emulsion PCR [1201]. This produces beads conjugated to dsDNA comprising multiple different regions. In some cases the dsDNA [1240] comprises seven regions: a first adaptor region [1241], a promoter region [1242], a partial sequencing adaptor region [1243], a randomer cell tag region [1244], a universal sequence [1245], a restriction site [1246], and a second adaptor region [1247]. The first and second adaptor regions can be 454 A and 454 B adaptors, respectively. The promoter region can be a T7 promoter. The partial sequencing adaptor region can be a partial Illumina adaptor. The universal sequence can end with GGG. A restriction digest may then be conducted. In some cases, the restriction digest is a blunt end restriction digest [1202]. This modifies the dsDNA bound to the bead surface. In some cases the surface-bound DNA then comprises five regions: an adaptor region [1251], a promoter region [1252], a partial sequencing adaptor region [1253], a randomer cell tag region [1254], and a universal sequence [1255]. The adaptor region can be a 454 A adaptor. The promoter region can be a T7 promoter. The partial sequencing adaptor region can be a partial Illumina adaptor. The universal sequence can end with GGG. In vitro transcription is then conducted [1203], producing RNA oligonucleotide products [1260]. In some cases, the oligonucleotide products comprise three regions: a partial sequencing adaptor region [1263], a randomer cell tag [1264], and a universal sequence [1265]. The partial sequencing adaptor region can be a partial Illumina adapter. The second universal sequence can end with GGG.

FIG. 13 depicts another method of producing beads conjugated to dsDNA comprising a cell tag. In some cases the dsDNA [1311] bound to the bead [1310] comprises one region: an adaptor region [1311]. The adaptor region can be a 454 A adaptor. In some cases, the oligonucleotides [1320] in solution comprise five regions: an adaptor region [1321], a promoter region [1322], a partial sequencing adaptor region [1323], a randomer cell tag region [1324], and a universal sequence [1325]. The adaptor region can be a 454 A adaptor. The promoter region can be a T7 promoter. The partial sequencing adaptor region can be a partial Illumina adaptor. The universal sequence can end with GGG. In some cases, the primers in the solution [1335] comprise a universal sequence. The universal sequence can end with GGG. PCR is then conducted. In some cases, the PCR is emulsion PCR [1301]. The PCR reaction may result in bead conjugated to dsDNA comprising a cell tag and the indicated regions. In some cases the dsDNA [1340] comprises five regions: an adaptor region [1341], a promoter region [1342], a partial sequencing adaptor region [1343], a randomer cell tag region [1344], and a universal sequence [1345]. The adaptor region can be a 454 A adaptor. The promoter region can be a T7 promoter. The partial sequencing adaptor region can be a partial Illumina adaptor. The universal sequence can end with GGG. The first and second adaptor regions can be 454 A and 454 B adaptors, respectively. The promoter region can be a T7 promoter. The partial sequencing adaptor region can be a partial Illumina adaptor. The universal sequence can end with GGG. In vitro transcription is then conducted [1302], producing RNA oligonucleotide products [1350]. In some cases, the oligonucleotide products comprise three regions: a partial sequencing adaptor region [1353], a randomer cell tag [1354], and a universal sequence [1355]. The partial sequencing adaptor region can be a partial Illumina adapter. The second universal sequence can end with GGG.

Example 2—High Throughput Single Cell RNA Tagging by Merging of Droplets

A first solution, containing a population of cells, is partitioned into a first series of droplets by a droplet generator (e.g. BioRad QX200 system, Dolomite Microfluidics systems, Micronit Microfluidics systems, water-in-oil microfluidic T-junction devices) using Poisson distribution statistics to ensure that each droplet contains 0 or 1 cells. A second solution, containing a population of beads from Example 1 as well as oligonucleotides, primers, enzymes, and other necessary reagents for PCR, in vitro transcription, and reverse transcription, is partitioned into a second series of droplets by a droplet generator using Poisson distribution statistics to ensure that each droplet contains 0 or 1 beads. Droplet generator systems are capable of generating tens of thousands of droplets, enabling high-throughput screening. Double-stranded DNA (dsDNA) is bound to the surface of the beads. The dsDNA may contain: (a) sequences for a promoter region (e.g., T7 promoter); (b) a sequencer adaptor region; (c) a randomer molecule tag region which is different for each DNA molecule on the bead; (d) a first universal sequence; (e) a second universal sequence ending with GGG; and/or (f) a randomer cell tag region which is the same for each DNA molecule on a particular bead. Droplets from the first series which contain 0 cells are discarded, and droplets from the second series which contain 0 beads are discarded. Discarding of droplets is performed by a droplet sorter (e.g. fluorescence-activated droplet sorting (FADS) system).

A set of reactions are carried out within the first series of droplets and the second series of droplets. The reactions may be conducted in parallel or sequentially. The first series of droplets is subjected to conditions which lyse the cells within the droplets. Reverse transcription is then conducted on the mRNA from the cell lysate, thereby producing cDNA transcripts within the droplets. The reverse transcription reaction may be performed using primers that are specific for certain regions of the mRNA; in some cases, the primers are general primers (e.g. as in QIAGEN QuantiTect Reverse Transcription Kit). (In some cases, the reverse transcription reaction occurs after a droplet from the first series of droplets is merged with another droplet, such as a droplet from the second series of droplets.) In the second series of droplets, the DNA molecules on the beads serve as templates for in vitro transcription reactions within the droplets which produce a set of mRNA transcripts from the DNA molecules encoding cellular and/or molecular tags. (In some cases, the in vitro transcription reaction occurs after a droplet from the second series is merged with another droplet, such as a droplet from the first series of droplets.)

Each single cell lysate-containing droplet is merged with one single bead transcript-containing droplet using a microfluidic device (e.g., RainDance RainDrop system, Dolomite Microfluidics systems). The GGG sequence in the mRNA tags then hybridizes to the overhanging CCC within the cellular cDNA. A reverse transcription reaction is then performed using the RNA tag as a template for an extension reaction. (e.g. FIG. 3, step [303]), resulting in tagging of the cellular cDNA with a DNA molecular/cellular tag.

Universal PCR is conducted on the pooled cDNA to further amplify the products. Target-specific PCR is conducted to target the genes of interest for sequencing. These target specific PCR products contain the genes of interest as well as the randomer cell and molecule tags and the sequencing adaptor sequence. These target specific PCR products are sequenced (e.g. by Illumina sequencing, where the sequencing adaptor sequences in the PCR products are Illumina tags). The information in the cell tag sequences is used to correlate specific genetic information with specific cells. The information in the molecule tag sequences is used to correlate specific genetic information with specific original mRNA molecules from each cell, allowing for correction of PCR amplification bias.

Example 3—Measurement of the Immune Response

A functional immune response (IR) assay is performed on sample comprising a population of cells, each expressing different immune receptors and comprising genes encoding immunoglobulin heavy chain or TCR beta, immunoglobulin light chain or TCR alpha, and various other genes related to immune function. The sample is partitioned into partitions containing immune cells. Each partitioned cell population may be lysed and its RNA may be extracted. The RNA from cell population is tagged with a molecule tag where the molecule tag is different for each molecule from a given cell by performing reverse transcription using an immunoglobulin or T-cell receptor and a template switching oligonucleotide comprising a molecule tag, and adapter sequence.

After tagging, the tagged RNA or its corresponding cDNA is pooled, amplified, and sequenced for analysis. The information from the expressed RNA may be used in combination with information from the molecule tags to determine the diversity of the immune response of a subject

Example 4—Measurement of the Immune Response at the Single Cell Level

A functional immune response (IR) assay is performed on sample comprising a population of single cells ([611], [612], [613], [614], [615], [616], [617], [618], [619]), each expressing different immune receptors [603] and comprising genes encoding immunoglobulin heavy chain or TCR beta [601], immunoglobulin light chain or TCR alpha [602], and various other genes related to immune function [604]. The sample is partitioned into partitions containing no more than one cell. Each partitioned cell may be lysed and its RNA may be extracted [605]. The RNA from each single cell is tagged with a cell tag where the cell tag is the same for every molecule from one given cell, and is also tagged with a molecule tag where the molecule tag is different for each molecule from a given cell.

After tagging, the tagged RNA or its corresponding cDNA ([621], [622], [623], [624], [625], [626], [627], [628], [629]) is pooled, amplified, and sequenced [606] for analysis. The information from the expressed RNA may be used in combination with information from the cell and/or molecule tags to categorize the cells into groups based on their immune function gene expression ([631], [632], [633], [634], [635], [636], [637], [638], [639]). This captures the information traditionally provided by multiple assays. For example, this immune response assay can provide information akin to that from an immune repertoire assay [710], flow cytometry [720], and a gene expression assay [730].

Example 5—Single Cell cDNA Tagging

FIG. 14 depicts incorporation of a cell tag onto the 5′ end of a cDNA produced from cellular mRNA (corresponding to the 3′ end of the mRNA) as well as incorporation of a molecular tag onto the 3′ end of the cDNA (corresponding to the 5′ end of the mRNA). The molecular tag is incorporated using RACE in this case and the cellular tag is added using a cDNA synthesis primer. The same procedure can be used to incorporate the cell tag onto the 3′ end of the cDNA and the molecular tag onto the 5′ end of the cDNA.

One cell or lysate from one cell is placed in a compartment. The cell may then be lysed in the compartment if it was not previously lysed. The lysate contains mRNA molecules [1410] containing an mRNA sequence and poly(A) tail [1414]. The compartment also contains oligodT molecules [1420]. In some cases, the oligodT molecules comprise five regions: an NVT region [1425], a TCA region [1426], a first universal sequence [1427], a randomer cell tag region [1428] where all oligodT molecules in a compartment have the same cell tag sequence, and a second universal sequence [1429]. The compartment also contains oligonucleotides [1430] which may comprise three regions: a universal sequence [1431], a randomer molecule tag [1432] where all oligonucleotides in the compartment have a different molecule tag sequence, and an overhang region [1433]. A reverse transcription reaction may be conducted using the oligodT region [1425] to primer the polyA tail present on the mRNA [1410]. The template switching oligonucleotide [1430] comprising the molecule tag may then be used as a template for the 5′ RACE assay to add the molecule tag to the cDNA. This results in cDNA molecule products [1440]. In some cases, the cDNA molecules contain nine regions: a first universal sequence [1441], a randomer molecule tag [1432] where all original oligonucleotides in the compartment have a different molecule tag sequence, an overhang region [1433], an mRNA information sequence [1444] containing the information from the original mRNA molecules, an NVT region [1435], a TCA region [1436], a second universal sequence [1437], a randomer cell tag region [1438] where all cDNA molecules in a compartment have the same cell tag sequence, and a third universal sequence [1439].

The labeled cDNA products in a compartment are then pooled with labeled cDNA products from other compartments and amplified with universal PCR. The amplified cDNA products [1510] may comprise nine regions: a first universal sequence [1511], a randomer molecule tag [1512] where all original oligonucleotides from a given compartment have a different molecule tag sequence, an overhang region [1513], an mRNA information sequence [1514] containing the information from the original mRNA molecules, an NVT region [1515], a TCA region [1516], a second universal sequence [1517], a randomer cell tag region [1518] where all cDNA molecules from a given compartment (and therefore the products from those molecules) have the same cell tag sequence, and a third universal sequence [1519]. These PCR products are melted and the single stranded DNA (ssDNA) is circularized with enzymes. The circularization product cDNA molecules may comprise nine regions: a first universal sequence [1521], a randomer molecule tag [1522] where all original oligonucleotides from a given compartment have a different molecule tag sequence, an overhang region [1523], an mRNA information sequence [1524] containing the information from the original mRNA molecules, an NVT region [1525], a TCA region [1526], a second universal sequence [1527], a randomer cell tag region [1528] where all cDNA molecules from a given compartment (and therefore the products from those molecules) have the same cell tag sequence, and a third universal sequence [1529]. The circularization allows reordering of the regions to group together regions of interest for sequencing. Target specific PCR is conducted to amplify regions containing genes of interest as well as cell and molecular tags, which may then be sequenced and the information processed. The target specific PCR products may comprise seven regions: a first universal sequence [1531], a randomer molecule tag [1532] where all original oligonucleotides from a given compartment have a different molecule tag sequence, an overhang region [1533], an mRNA information sequence [1534] containing the information from the original mRNA molecules, a second universal sequence [1537], a randomer cell tag region [1538] where all cDNA molecules from a given compartment (and therefore the products from those molecules) have the same cell tag sequence, and a third universal sequence [1539].

Example 6—Single Cell cDNA Tagging Using a Microwell Array

The single cell cDNA tagging described in Example 4 may be conducted using the wells on a first microwell array plate as the compartments [1610], as depicted in FIG. 16. The oligodT molecules may be coupled to beads ([1611], [1612], [1613], [1614], [1615], [1616]) which are deposited with one bead per well. The microwell array plate may be designed to permit only one bead per well. For example, the wells may be no more than 100 micron in diameter. Master cDNA mix comprises the necessary reagents and enzymes, including the oligonucleotides with randomer molecule tags, for the generation of the labeled cDNA products described in Example 4. This master cDNA mix is also deposited in the wells. Random template-switching oligonucleotides may also be deposited into the wells. The wells may also optionally contain a lysis reagent. A second microwell array plate [1620] is loaded with cells ([1621], [1622], [1623], [1624], [1625], [1626]), with no more than one cell per well. The microwell array plate may be designed to permit only one cell per well. For example, the wells may be no more than 30 micron in diameter. The cells may be in PBS solution. The two microwell array plates are aligned to bring individual wells on the first plate in contact with individual wells on the second plate [1601]. The wells may be placed in an incubator at 42° C. Contact between the cell in one well and the cDNA mix in one other well causes lysis of the cell [1602]. After lysis, the beads may be collected and universal PCR may be conducted on the beads. In some instances, the compartments formed by the matched microwells are used as the compartments for the process described in Example 4.

Example 7—Single Cell cDNA Tagging Using a Microarray

The single cell cDNA tagging described in Example 4 may be conducted using the spots on a microarray, partitioned from each other, as the compartments [1640]. The oligodT molecules may be bound to individual spots on the array ([1641], [1642], [1643], [1644], [1645], [1646]) which are then partitioned from each other to form compartments. In some cases, the partitioning is achieved by aligning a PDMS membrane with wall or chamber features onto the array surface. Master cDNA mix comprises the necessary reagents and enzymes, including the oligonucleotides with randomer molecule tags, for the generation of the labeled cDNA products described in Example 4. This master cDNA mix is deposited hi the chambers. A microwell array plate [1650] is loaded with cells ([1651], [1652], [1653], [1654], [1655], [1656]), with no more than one cell per well. The cells may be in PBS solution. The microwell array plate and the microarray are aligned to bring individual wells on the microwell array plate in contact with individual chambers on the second plate [1601]. Contact between the cell in one well and the cDNA mix in one other well causes lysis of the cell [1603]. The compartments formed by the matched microwells and microarray chambers are then used as the compartments for the process described in Example 4. 

What is claimed is:
 1. A method of tagging target oligonucleotides comprising: a. partitioning DNA into a plurality of compartments; b. performing an in vitro transcription reaction on the DNA within the compartments, thereby obtaining compartments comprising RNA; c. merging the interior of the compartments comprising RNA with the interior of a set of compartments comprising target oligonucleotides; d. hybridizing the RNA to the target oligonucleotides; and e. performing a reaction to attach a sequence corresponding to the RNA to the target oligonucleotides.
 2. The method of claim 1, wherein the DNA is double stranded.
 3. The method of claim 1, wherein the compartments are droplets within an oil-and-water emulsion.
 4. The method of claim 1, wherein the target oligonucleotides comprise at least one target oligonucleotide comprising a cell tag and a molecule tag.
 5. The method of claim 1, wherein the target oligonucleotides are DNA.
 6. The method of claim 5, further comprising, before the merging of step (c), partitioning the target oligonucleotides into the set of compartments.
 7. The method of claim 5, further comprising, before the merging of step (c), partitioning a set of cells into the set of compartments and lysing the cells in order to liberate cellular oligonucleotides.
 8. The method of claim 7, wherein the cellular oligonucleotides are the target oligonucleotides.
 9. The method of claim 7, wherein the cellular oligonucleotides are cellular mRNA.
 10. The method of claim 9, further comprising conducting a reverse transcription on the cellular mRNA to generate cellular cDNA.
 11. The method of claim 10, wherein the target oligonucleotides are cellular cDNA.
 12. The method of claim 1, wherein the reaction is a Rapid Amplification of cDNA Ends (RACE) reaction.
 13. The method of claim 10, wherein the reverse transcription reaction is performed with a primer specific for a region of the genome.
 14. The method of claim 1, wherein the DNA is conjugated to a solid support.
 15. The method of claim 14, wherein the solid support is a bead.
 16. The method of claim 10, wherein the reverse transcription reaction is conducted within the set of compartments prior to the merging step.
 17. The method of claim 10, wherein the reverse transcription reaction is conducted within the merged compartments.
 18. The method of claim 13, wherein the region of the genome is an immunoglobulin gene or a T-cell receptor gene.
 19. A method comprising: (a) providing a plurality of beads comprising a plurality DNA oligonucleotides; (b) providing a plurality DNA oligonucleotides comprising a primer sequence, a universal sequence an adapter sequence and a cellular tag; (c) merging the beads of step (a) and the oligonucleotides of step (b) into a plurality of compartments such that each compartment comprises a single bead and a single oligonucleotide; (d) performing an amplification reaction on the oligonucleotides within the compartments, thereby obtaining a plurality of DNA oligonucleotides comprising the primer sequence, the universal sequence, the adapter sequence and the cellular tag; (e) performing an in vitro transcription reaction on the DNA within the compartments, thereby obtaining compartments comprising RNA the primer sequence, the universal sequence and cellular tag; merging the interior of the compartments comprising RNA with the interior of a set of compartments comprising target oligonucleotides; (g) hybridizing the RNA to the target oligonucleotides; and (h) performing a reaction to attach a sequence corresponding to the RNA to the target oligonucleotides.
 20. The method of claim 19, wherein the plurality of oligonucleotides on the bead comprise a molecule tag.
 21. A method of tagging target oligonucleotides comprising: (a) isolating a plurality of mRNA from a biological sample comprising a plurality of cell types; and (b) performing reverse transcription of the mRNA using a primer specific for the target oligonucleotide and a template switching oligonucleotide comprising a molecule tag, a universal sequence, and an adapter sequence to produce tagged target cDNA
 22. The method of claim 21, wherein the target cDNA is tagged at the 3′end.
 23. The method of claim 21, wherein the target oligonucleotide is an immunoglobulin or T-cell receptor.
 24. The method of claim 21, wherein the adapter sequence is specific to a sequencing platform.
 25. The method of claim 21, wherein the molecule tag is an oligomer.
 26. The method of claim 25, wherein the oligomer is a randomer.
 27. The method of claim 26, wherein the randomer is at least a 9mer.
 28. The method of claim 21, further comprising amplifying the target cDNA using the universal sequence and a primer specific for the target oligonucleotide.
 29. The method of claim 28, further comprising sequencing the amplified cDNA.
 30. A method of determining the immune repertoire in a subject comprising; (a) isolating a plurality of mRNA from a biological sample comprising a plurality of cell types; and (b) performing reverse transcription of the mRNA using a immunoglobulin or T-cell receptor specific primers and a template switching oligonucleotide comprising a molecule tag, a universal sequence and an adapter sequence to produce molecule tagged immunoglobulin or T-cell receptor cDNA; (c) amplifying the cDNA using the universal sequence and a primer specific for the target oligonucleotide; (d) sequencing the cDNA to produce a plurality of sequencing reads; (e) grouping the sequence reads with the same molecule tag and clustering the sequences within the same group; and building a consensus sequence for each cluster to produce a collection of consensus sequences wherein the consensus sequence is used to determine the diversity of the immune repertoire.
 31. The method of claim 30, wherein the target cDNA is tagged at the 3′ end. 